
function xOut=times(x1,x2)
  n1 = size(x1);
  n2 = size(x2);
  x = getval(x1);
  y = getval(x2);
  dfdx1 = y;
  dfdx2 = x;
  dfdx1 = dfdx1(:);
  dfdx2 = dfdx2(:);
  if(notisa(x1,'deriv1s'))
    np = nindep(x2);
  else
    np = nindep(x1);
  end
  if(notisa(x1,'deriv1s'))
    xOut.v = times(x1,x2.v);
    xOut.d = elelmult_splastindx(dfdx2,x2.d);
  elseif(notisa(x2,'deriv1s'))
    xOut.v = times(x1.v,x2);
    xOut.d = elelmult_splastindx(dfdx1,x1.d);
  else  % both deriv1
    n1d = size(x1.d);
    np = n1d(end);
    xOut.v = times(x1.v,x2.v);
    xOut.d = elelmult_splastindx(dfdx1,x1.d) + elelmult_splastindx(dfdx2,x2.d);
  end
  xOut=class(xOut,'deriv1s');
